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Using scanning tunneling spectroscopy (STS), we address the problem of the superconductor- 
insulator phase transition (SIT) in homogeneously disordered ultrathin (2-15 nm) films of NbN. 
Samples thicker than 8 nm, for which the loffe-Regel parameter kpl > 5.6, manifest a conventional 
superconductivity : A spatially homogeneous BCS-like gap, vanishing at the critical temperature, 
and a vortex lattice in magnetic field. Upon thickness reduction, however, while kpl lowers, the 
STS revealed striking deviations from the BCS scenario, among which a progressive decrease of 
the coherence peak height and spatial inhomogeneities. The thinnest film (2.16 nm), while not 
being exactly at the SIT ( Tc ~ OATc-buik), showed astonishingly vanishing coherence peaks and 
the absence of vortices. In the quasi- 2D limit, such clear signatures of the loss of long-range phase 
coherence strongly suggest that, at the SIT the superconductivity is destroyed by phase fluctuations. 
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Despite numerous theoretical and experimental works 
[l| , the understanding of the phase transition from a su- 
perconducting (SC) to an insulating state in ultrathin 
films remains a very challenging problem. The struc- 
tural properties of the films are known to play a key-role 
in determining the low-temperature electronic properties. 
Disordered SC thin films undergoing a superconductor- 
insulator transition (SIT) can be divided into two groups 
Q : i) Granular systems, consisting of coupled SC grains, 
described at low temperature as an array of Josephson 
junctions [3|; ii) Systems qualified as "homogeneous" 
where both amplitude and phase variations of the or- 
der parameter are believed to play an important role Q . 
In contrast, in perfectly crystallized systems the SIT was 
not observed down to a single atomic monolayer [5|. 

A way to get a new insight in the microscopic pro- 
cesses occurring at the SIT is to study the local properties 
of the SC condensate evolving in a disordered potential. 
Recently two scanning tunneling microscopy and spec- 
troscopy (STM/STS) experiments addressed this ques- 
tion bringing new interesting results and leading to var- 
ious successive interpretations [sl, Q • It is worth noting 
that the SIT problem has also important application is- 
sues, specifically in the field of single photon detectors 
m in which ultrathin SC films (usually NbN or NbTiN) 
are used. While the microscopic picture of the photon- 
to-condensate excitation conversion does not imply any 
thickness dependence, it has recently been discovered 
that the maximum detection efficiency occurs for film 
thicknesses close to the SIT [9|. Thus, the detailed study 
at the nanometer scale of such SC films would be of help 
in the understanding of the underlying microscopic pro- 
cesses of the photon absorption by real SC devices. 



In this Letter, we address the SIT problem locally in a 
direct STS experiment by varying the thickness of ultra- 
thin NbN films from 15 down to 2.16 nm, the correspond- 
ing Tc varying from 15.0 to 6.7 K, while the loffe-Regel 
parameter kpl decreases from 5.7 to below 2.6. Our STS 
data reveal profound changes in the local behavior of the 
SC films as the SIT is approached. The amplitude of 
the coherence peaks diminishes to almost vanish when 
Tc ^ OATc-buik^ which is not in the very vicinity of 
the insulating transition. This feature is accompanied by 
a significant broadening of the quasiparticle density of 
states (DOS) and by the development of a correlated elec- 
tronic background evidenced by a characteristic V-shape, 
extending to voltages 5-10 mV significantly larger than 
A/e. In addition, a shallow pseudogap regime is observed 
in a very narrow temperature window just above Tc- 
Interestingly, our results on 2D-superconductors present 
many similarities with the results reported for 3D-NbN 
films of comparable Tc and kpl [1Q| but are in strik- 
ing contrast with the recently reported large pseudogap 
regimes observed in TiN and InO [sl, • We ascribe these 
differences to a more homogeneous crystalline structure 
of our NbN films and to film conductances farther away 
from the critical one. Finally, while the vortex lattice 
is observed in thicker films subject to magnetic field, the 
STS contrast decreases when reducing the film thickness; 
strikingly, the 2.16 nm film with Tc ^ OATc-buik does 
not show any evidence of the vortex state, emphasizing 
thereby the loss of long-range SC phase coherence. Our 
results on 2D-homogeneously disordered superconductors 
therefore strongly support that phase fluctuations drive 
the insulating transition and have already profound ef- 
fects before the SIT. 
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FIG. 1: (Color online) a) Square resistivity vs temperature 
dependence of the studied samples; the dashed line indicate 
h/4e^ the quantum resistance for pairs; Tc decreases with 
film thickness (T^i^^"" = 15.0 K (not shown), T^S^^ = 14.5 K, 

^4nm ^ -^3 3 ^ ^2.33nm ^ g ^ ^2.16nm ^ ^y^ 

defined Tc by the temperature reached when the resistivity 
equals 10% of its 40 K value, b) Typical topographic STM 
image of the surface (15 nm-thick film, T = 4.2 K, image size 
l/xm X l/xm); c) The evolution of the tunneling conductance 
spectra (thin black lines) with temperature follows a conven- 
tional BCS behavior in the 15 nm sample: a SC gap around 
zero-bias, A(0) = 2.85 meV, which vanishes at the critical 
temperature of the film, with = 4.4; the bottom thick 

blue line is a BCS fit at 2.3 K with no additional broadening 
parameter, d) dl/dV conductance map at = in magnetic 
field (B — 1 T), revealing a disordered vortex lattice in the 
15 nm film. 



In order to achieve as realistic conditions as possible, 
ultrathin NbN films used in real single photon detectors 
[sl were probed; for every thickness the films were grown 
using identical growth conditions on crystalline sapphire 
substrates by DC reactive magnetron sputtering of a pure 
Nb target in an Ar+N2 gas mixture. Their structural, op- 
tical, electronic and SC properties were extensively stud- 
ied and are thus well established [ll|; 12 nm thick films 
have a Tq about 15 K and kpl ^ 5.7; these parameters 
hardly change for 6 nm thick films for which Tc is about 
14 K and kpl ^ 5.6; reducing further the thickness by 
two, 3.2 nm, Tc is reduced to 10.7 K and kpl ^ 2.6; the 
structural characterization shows that such NbN films are 
covered by a 0.5-1 nm thin natural passivation layer [Tlj . 
Estimation of the Ginzburg-Landau coherence lengths re- 
veals values between 4 and 6 nm; thus films having thick- 
nesses below about 5 nm can be considered as true 2D 
SC, the passivation layer having to be considered as be- 
longing to the system itself. For thicker films, we see from 



the Tc values very close to bulk Tc ^ 15.5 K, that the 
passivation layer induces a negligible inverse proximity 
effect on SC. 

Rsq{T), the square resistivity vs temperature depen- 
dence of the studied films, is presented in Fig. la. It shows 
a typical SIT behavior of non-granular films : Rsq{T) 
keeps a very smooth dependence on approaching Tc; Tc 
monotonically lowers as the film thickness decreases. Be- 
ing weak for thicker films (8 nm and more), this effect 
becomes significant for thinner films, reflecting the rapid 
suppression of Tc on approaching the SIT. Among the 
studied samples, the 2.16 nm one is the closest to the SIT, 
"sitting" on its SC side, and having Tc ^ OATc-buik] yet 
thinner films (2.0 nm) show an insulating behavior at low 
temperature. 

All the studied films were grown ex situ. No fur- 
ther treatment was performed prior to STM/STS exper- 
iment. Constant current STM images revealed a very 
smooth film surface with ^ 0.1 nm roughness, which 
is 10 times less than the corrugation of patterned films 
The STS was performed at variable temperatures 
(2.5-20 K) using mechanically cut Pt/Ir tips. The local 
tunneling conductance spectra dI/dV{V^x^y) were ob- 
tained by numerical derivative of raw I{V) data. They 
are directly linked to the local quasiparticle excitation 
spectrum, A^^ {E^ x, ?/), through the relation dI{V)/dV ex 
- / dENs {E + eV,x,y)df{E)/dE, where f{E) is the 
Fermi-Dirac function. 

The thicker film, 15 nm, shows conventional and 
spatially homogeneous SC conductance spectra, with a 
non-broadened BCS-like shape, characterized by a gap 
A(0) = 2.85 meV (Fig.l); the gap was observed to van- 
ish at the film's Tc extracted from Rsq(T) measurement. 
Upon application of a magnetic field, a well-known disor- 
dered vortex lattice was revealed, attesting for local dis- 
order existing in this extreme type II SC (Fig. Id). Thus, 
our STS data on the 15 nm sample are in goo d agree ment 
with the conclusions of previous reports [3, EMI. 

Figure 2 summarizes our STS results obtained on thin- 
ner films (8, 4, 2.33 nm). The data are organized in three 
columns by thickness. From top to bottom, each column 
shows: i) a normalized conductance map for a bias volt- 
age chosen to picture out possible gap inhomogeneities, 
ii) selected spectra representative of the maximal varia- 
tions observed among the local dl/dViV) spectra (red 
and blue curves) and typical spectrum in magnetic field 
at vortex centers (black curve), iii) an image of the vortex 
lattice and iv) the temperature evolution of local spectra. 

The 8 nm dl/dV map reveals only tiny spatial inho- 
mogeneities of the SC gap, except in regions where struc- 
tural defects occur, such as the red patch in the image. 
The disordered vortex lattice is routinely observed. The 
temperature evolution of the tunneling spectra remains 
conventional, although a slight spectral broadening was 
detected and accounted for with a very small Dynes pair- 
breaking parameter F ^ 0.01 meV. In overall, the 8 nm 
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FIG. 2: (Color online) From left to right column: tunnel- 
ing characteristics of 8, 4 and 2.33 nm thick films. From 
top to bottom: the first row shows normalized STS conduc- 
tance maps representing the color-coded spatial variations of 
the dl/dV signal measured at 2.3 K at the gap onset (for 
V = 1.7 mV, V = 1.7 mV and V = 1.1 mV); second row: 
normalized individual tunneling conductance spectra repre- 
sentative of the spatial variations observed among the local 
conductance spectra (red and blue curves) at 2.3 K; the black 
curves show typical spectra measured at 4.0 K in magnetic 
field {B = 1 T, B = 1 T and B = 2 T) at the vortex centers; 
third row: dl/dV conductance maps at V = showing im- 
ages of the vortex lattices at 4.0 K (for B = 1T,B = 1T and 
5 = 2 T); fourth row: temperature evolution of the normal- 
ized dI/dV{V) tunneling spectra (vertically shifted for clar- 
ity); the temperature is indicated below the corresponding 
spectra. Red (black) curves correspond to the SC (normal) 
state. The bottom thick blue lines are BCS dl/dV calcula- 
tions at the indicated temperature; an additional broadening 
parameter was needed to fit the data (see text). 



NbN film behaves similarly to the 15 nm one. 

The situation changes in the 4 nm sample. While the 
vortex lattice and the temperature evolution of the tun- 
neling spectra do not manifest significant differences with 
respect to thicker samples, small spatial inhomogeneities 
are revealed in the tunneling conductance map, which 
specifically affect the gap width and the height of the 



coherence peaks (see the red and blue curves below the 
conductance map). These inhomogeneities are charac- 
terized by two spatial length scales. The smaller one is 
of the order of a few nm, thus close to the coherence 
length scale. The larger one is of several tens of nm; fur- 
ther analysis showed that the latter originates from slight 
variations of the film thickness due to underlying atomic 
steps of the substrate (see also Fig. 3 and section 3 of the 
supplemental material). Nevertheless, the dl/dV spec- 
tra at T <C Tc remain reasonably BCS-like, the fit again 
requiring F 0.01 meV. 

In the 2.33 nm sample, the tunneling spectra undergo 
more and more pronounced changes. The conductance 
maps reveal inhomogeneities comparable to the 4 nm 
case, characterized by two length scales. "Full" gap in 
dl/dV exists everywhere, characterized by zero conduc- 
tance and coherence peaks, with however their amplitude 
much reduced with respect to the case of thicker films; 
their height is observed to vary spatially at the nanometer 
scale (see the blue and red curves), similarly to the effects 
discovered in cuprates [15|. The vortex lattice is hardly 
observable; Inside the cores the spectra present a large 
dip with no coherence peaks. Moreover, the normal state 
spectra reveal a V-shaped background with a clear mini- 
mum at zero-bias, extended on a larger energy scale than 
the SC gap. A shallow pseudogap (a dip opening above 
Tc) is also observed within a 1 K range above Tc, whose 
energy scale is comparable to the SC gap. The spectra 
are poorly described by a BCS expression at T <C Tc and 
require F ^ 0.05 meV (see fit). 

Striking changes are revealed in the tunneling charac- 
teristics of the 2.16 nm sample (see Fig. 3). The SC gap 
develops on a much deeper V-shaped background, exist- 
ing also above Tc (Fig. 3e) thus characteristic of the nor- 
mal state, the energy scale of which is of 5-10 meV, sim- 
ilar to the one observed for the 2.33 nm sample. In addi- 
tion, a very shallow pseudogap is present; it appears more 
clearly in Fig. 3f-g where Rsq{T) and dl/dvlv = 0,T) 
are presented simultaneously for 2.16-8 nm films. The 
difference between the two behaviors is straightforward 
: In thicker samples, the SC gap opens at Tc on a flat 
spectral background (usual metallic state), while in the 
thinner 2.16 and 2.33 nm samples, the DOS at Ep is 
already reduced in the normal state. Additionally, the 
opening of the SC gap occurs at Tc, as it is unambigu- 
ously detected as a kink in the dI/dV{V = 0,T) curves 
(pointed by the arrows), in contrast to the smoother evo- 
lution in high-Tc cuprates (see [16] and ref. therein). 

Furthermore, the spectra are not BCS-like anymore; 
using nevertheless a broadened BCS DOS to estimate 
the gap value, yields F ^ 0.10 meV. The amplitude of the 
quasiparticle peaks is further damped with respect to the 
2.33 nm case, so that in some areas the peaks disappear 
completely. The conductance map in Fig. 3a allows again 
visualizing two-scale spatial variations, as observed in 4 
nm and 2.33 nm samples, with spectral inhomogeneities 
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FIG. 3: (Color online) Tunneling characteristics of the 
2.16 nm thick film at 2.3 K. a) zero-bias STS conductance 
maps at B = T, b) at B = 3 T. c) Typical variations 
among the corresponding local spectra for B = T and d) 
for B = 3 T; e) Evolution of the dlldV(y) spectra with 
temperature. The bottom thick blue line is a BCS calcula- 
tion at 2.3 K with an additional broadening parameter (see 
text), f) zoom on i?(T) of the 2.16-8 nm films compared to 
their g): temperature dependence of the zero-bias conduc- 
tance dliy ^^)ldV. 



quite comparable to these latter cases, i) The zero-bias 
conductance and the height of quasiparticle peaks vary 
slightly at the nm scale, ii) The large scale variations 
appear strikingly to be spatially ordered in parallel bands 
in some areas of the sample (Fig. 3a-b; see section 3 of 
the supplemental material for more details); 

Another astonishing feature is observed : the spec- 
tra are dramatically unsensitive to the magnetic field 
(Fig. 3b-d); up to B=3 T, no vortices could be observed 
(Fig. 3b). However, while bulk Rc2 is estimated to be 

20 T, lic2 decreases with film thickness, but remains 
larger than 10 T for all our studied films and is thus 
much higher than the fields used here. Vortices are a 
mesoscopic manifestation of the quantum coherence of 



TABLE I: Critical temperature Tc (extracted from transport 
measurements), energy gap value Abcs (extracted from BCS 
fits of the tunnel conductance spectra) and 2Abcs /ksTc ra- 
tio as a function of film thickness d. Although the tunnel con- 
ductance spectra for the thinnest films at the lowest studied 
temperature (2.3 K) deviate from BCS behavior, we estimate 
the energy gap from a broadened BCS simulation. 
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the SC condensate, where the phase of the order param- 
eter is spatially well-defined and makes 27r turns around 
each vortex core. Thus, we interpret this striking result 
as a break-up of the long-range phase coherence in this 
fragile SC condensate, suggesting a lowered phase stiff- 
ness already at T = OATc-buik^ thus quite before the 
very SIT. 

A hallmark of our tunneling data on approaching the 
SIT is the development of a more and more pronounced 
V-shaped background when the film thickness is reduced. 
This suggest enhanced Coulombic effects and is con- 
sistent with the fact that our Tc vs Rsq data follow 
the Finkelstein law [l8[ (see figure and section 1 of the 
supplemental material). The observed A and Tc pro- 
gressive reduction (see Tab.|I]), with smoothly decreasing 
2Abcs /kTc ratio, is consistent with previous reports on 
homogeneous disordered thin films close to the SIT Q 
and with the behavior found for much thicker NbN films 
by changing the disorder [lo|. However, it is at variance 
with recent STS reports on TiN and InO films, both sug- 
gesting a huge increase of this ratio with increasing dis- 
order (61,12]. 

It is instructive to make a detailed comparison of 
our local results on NbN with the recently reported 
STM/STS experiments on TiN (3.6-5 nm), InO (15 and 
30 nm) and thicker (50 nm) NbN films. 

The most important difference with respect to TiN and 
InO samples, is that a very deep pseudogap above Tc was 
reported, where at Tc, almost all the DOS is depressed at 
Ef 0,0- Additionally, it is important to note that these 
systems were closer to the SIT {Tc < OATc-buik)- In out 
case, all Tc > OATc-buik and the largest depressed DOS 
at Ef, obtained for the closest films to the SIT, is of 
35 — 50%, but is mostly due to the V-shaped background 
and not to the pseudogap effect. On the other hand, 
there are also similarities: Observation of small gap in- 
homogeneities and suppression of the coherence peaks 
with increasing disorder while the conductance spectra 
deviate more and more from BCS behavior. 

The evolution of our tunneling spectra with varying 
thickness present many surprising similarities with recent 
studies on thicker NbN films, where the thickness was 
fixed and the SIT was approached by tuning the disorder 
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(from kpl 1.2 - 10.1) [lO|. In particular, thicker NbN 
films present a Tq versus kpl dependence quite compa- 
rable to ours, and lead to tunneling characteristics sim- 
ilar to what is reported here. In particular, the features 
reported for Tc = 6 K (/cf^ = 1.6) and Tc = 4.1 K 
{kpl = 2.2) are comparable with those presented for our 
2.16 nm film {Tc = 6.7, kpl < 2.6) in terms of: i) a large 
broadening of the spectra leading to gap filling, ii) almost 
vanishing of the coherence peaks, iii) the existence of a 
shallow pseudogap above Tc- 

For InO films, it was argued in [7] that the disappear- 
ance of the coherence peaks occurs at Tc. This is a no- 
ticeable difference with our results, where the decay of 
the coherence peaks amplitude is observed to be gradual 
as a function of thickness and Tc reduction. The co- 
herence peaks amplitude almost vanishes in the SC state 
(thus below Tc) when Tc ^ OATc-buik^ together with 
the concomitant disappearance of the vortex lattice, in- 
dicating the break-up of long range phase coherence of 
the SC condensate. Parts of the observed significant dif- 
ferences between NbN and TiN/InO films may be related 
to a more homogeneous crystalline structure of the NbN 
system, which is preserved and controlled upon thick- 
ness reduction (the growth conditions being identical for 
all films), leading to a more continuous behavior of the 
electronic properties toward the SIT. 

In conclusion, the present work shows that already at 
Tc ^ OATc-buik^ thus quite before the SIT, phase fluctu- 
ations of the SC order parameter lead to the almost van- 
ishing of the coherence peaks, together with the disap- 
pearence of the vortex lattice, while the SC state at zero 
magnetic field show only slight spatial spectral inhomo- 
geneities. Therefore, numerical simulations of a strongly 
disordered 5- wave superconductor, showing that ampli- 
tude fluctuations alone are not sufficient to drive the SIT 
[2q| are compatible with our results and with more recent 
numerical study [2l| . Further investigations very close to 
the SIT are necessary to clarify the role and the nature 
of the pseudogap and of the V-shaped background we 
observed in ultrathin homogeneous NbN films. 

The authors acknowledge fruitful discussions with B. 
Sacepe. The work is supported in part by DFG Center 
for Functional Nanostructures under sub-project A4.3. 
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